arXiv:1504.07884v3 [cond-mat.str-el] 31 Aug 2015 


Hyperspherical theory of the quantum Hall effect: the role of exceptional degeneracy 

K. M. Daily/ R. E. Wooten/ and Chris H. Greene^ 

^Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA 

(Dated: September 1, 2015) 

By separating the Schrodinger equation for N noninteracting spin-polarized fermions in two- 
dimensional hyperspherical coordinates, we demonstrate that fractional quantum Hall (FQH) states 
emerge naturally from degeneracy patterns of the antisymmetric free-particle eigenfunctions. In the 
presence of Coulomb interactions, the FQH states split off from a degenerate manifold and become 
observable as distinct quantized energy eigenstates with an energy gap. This alternative classification 
scheme is based on an approximate separability of the interacting Wfermion Schrodinger equation 
in the hyperradial coordinate, which sheds light on the emergence of Laughlin states as well as other 
FQH states. An approximate good collective quantum number, the grand angular momentum K 
from if-harmonic few-body theory, is shown to correlate with known FQH states at many filling 
factors observed experimentally. 

PACS numbers: 31.15.xj, 73.43.-f, 73.43.Cd 


I. INTRODUCTION 

One of the most striking aspects of nonrelativis- 
tic quantum mechanics in more than one dimension 
is the remarkable implication of high degeneracy or 
near-degeneracy. Textbook examples include the sp- 
hybridization of chemical bonds and the degenerate Stark 
effect of excited hydrogen atoms. In the degenerate Stark 
effect, for instance, even an infinitesimally small exter¬ 
nal electric field selects energy eigenstates that are linear 
combinations of a finite set of degenerate zero-field states, 
and these are the same eigenstates that can be obtained 
by separating the field-free Schrodinger equation for that 
system in parabolic coordinates. 

A major development presented in this article is that 
similar considerations apply to the fractional quantum 
Hall effect (FQHE) [H, I2|- In the fractional quantum 
Hall effect, a strongly-interacting two-dimensional elec¬ 
tron gas exhibits quantization in the presence of a strong, 
perpendicular magnetic field. In the noninteracting limit, 
the electrons fall into highly degenerate Landau levels, 
and their collective behavior depends on the filling fac¬ 
tor, the ratio of the number of electrons to the large, but 
finite, degeneracy of the lowest Landau level for a sample 
with finite area. The system exhibits quantization when 
the filling factor takes on integer or certain rational frac¬ 
tion values. We show that the high degeneracy of the 
noninteracting system produces dramatic implications. 

Moreover, we demonstrate that the N-electron 
Schrodinger equation is approximately separable in hy¬ 
perspherical coordinates. This approach to the problem 
shows that a characteristic property of the noninteract¬ 
ing system, which we denote the exceptional degeneracy, 
becomes unusually high for precisely those states that ap¬ 
pear at experimentally and theoretically observed FQHE 
filling factors. In other words, even though the FQHE 
is viewed fundamentally as the epitome of a strongly- 
correlated system of electrons, the occurrence or nonoc¬ 
currence of a FQHE filling factor is highly correlated 
with the pattern of exceptional degeneracies in the non¬ 


interacting electron system. The approximate separabil¬ 
ity demonstrated in hyperspherical coordinates for the 
few-body quantum Hall states makes definite predictions 
about a class of excitation frequencies that could be used 
to experimentally probe the system. 

Extensive progress in the theoretical understanding of 
the FQHE has been achieved through various approaches 
following the early intuitive development by Laughlin Q , 
notably the composite fermion (CF) picture developed 
by Jain [1,@, and work by Haldane 0, and Halperin Q. 
Theoretical treatments have tended to reside in one of 
two different categories, either postulating trial wave- 
functions as in Laughlin’s original approach, (e.g. 

0) or else performing large numerical diagonalizations 
for the maximum number of electrons that can still give 
a ma nag eable size computation, typically 8-20 particles 
(e.g. [il|). A more recently developed technique uses CF 
wave functions as a basis for numerical diagonalization to 
study systems with larger numbers of particles 

The approach developed here has some advantages 
complementary to previous methods. In contrast to tech¬ 
niques that use the single-particle representation (i.e. 
Slater determinant constructions), the approach treated 
here inherently uses collective coordinates. It also pro¬ 
vides a systematic expansion that can in principle de¬ 
scribe any states existing in the Hilbert space of a finite 
number of particles, while at the same time allowing us to 
see many key properties analytically or with small-scale 
diagonalizations. The adiabatic hyperspherical represen¬ 
tation capitalizes on an approximate separability, and its 
key element is a set of potential energy curves showing 
at a glance the relevant size and energy of different en¬ 
ergy eigenstates. While the hyperradial degree of free¬ 
dom is not separable for arbitrarily strong Coulomb inter¬ 
actions, our calculations demonstrate that approximate 
Born-Oppenheimer separability is an excellent approxi¬ 
mation in typical regimes of electron density and field 
strength for a typical material like GaAs. 

While the adiabatic hyperspherical representation di¬ 
ll^ has not been used extensively in condensed matter 
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theory, it has had extensive success in a wide range of 
few-body contexts. The literature in this field documents 
theoretical results that have been achieved in contexts as 
diverse as nuclear structure and reactivity uni¬ 

versal Efimov physics in cold atoms and molecules [2^ 
[ 23 , few-electron atoms and systems contain¬ 
ing positrons and electrons Efimov’s predic¬ 

tion of a universal binding mechanism for three 
particles at very large scattering lengths can itself be 
viewed as an application of the adiabatic hyperspheri- 
cal coordinate treatment in a problem where the method 
is exact, although Efimov did not himself express it in 
those terms. Hyperspherical coordinates have also been 
employed to describe some many-body phenomena such 
as the trapped atom Bose-Einstein condensate with ei¬ 
ther attractive or repulsive interactions (ssl . l3^ and the 
trapped degenerate Eermi gas in three dimensions includ¬ 
ing the BCS-BEC crossover problem. [H, 113 


Our initial presentation of the formulation begins by 
setting up the problem rigorously for N electrons con¬ 
fined to a plane with a transverse uniform magnetic field. 
The antisymmetric states for spin-polarized fermions 
are found directly using the technique developed in 
Ref. [13 • Next we show that the exact separability of the 
Schrodinger equation in hyperspherical coordinates for 
the case of noninteracting electrons still exhibits an ap¬ 
proximate separability even in the presence of Coulomb 
interactions. The treatment then demonstrates how in¬ 
teractions single out potential energy curves or channels 
of exceptional degeneracy, which correlate with known 
Laughlin states, composite fermion EQHE states, and 
suggest other states that deserve future theoretical and 
experimental investigation. This enables further predic¬ 
tions of a class of excitation frequencies that should be 
experimentally observable in a EQHE experiment. 


This paper is organized as follows. Sections HIl and 
m formulate the one-body and N-body relative Hamil¬ 
tonians in the symmetric gauge. Section IIVI defines the 
hyperspherical coordinates adopted in the present study 
and writes the unsymmetrized hyperspherical harmonics 
which serve as our primitive basis set. This basis set is 
then connected to the Landau level picture and suggests a 
definition for the hyperspherical filling factor. The effect 
of Coulomb interactions is then developed and treated 
within the adiabatic hyperspherical representation. Sec¬ 
tion El introduces the concept of exceptional degeneracy 
and computes this key quantity, which correlates with 
filling factors that are observable as EQH ground states. 
Section ED offers concluding remarks and comments on 
future directions. Einally, Appendix relates the states 
with different hyperspherical filling factors to several of 
the states that have previously been identified in the con¬ 
ventional composite fermion picture. 


II. SINGLE PARTICLE HAMILTONIAN, 


The noninteracting Hamiltonian H for a single electron 
in an external magnetic field is given by 

H=^(-ihV + eAf ( 1 ) 

ZrUe 


in SI units where nif, is the effective mass of the electron 
in the medium, e is the magnitude of the electron charge, 
h is Planck’s constant, and A is the vector potential. 
For two dimensional space, in Cartesian coordinates, the 
gradient is V = xdx + ydy. For a constant magnetic field 
of magnitude B oriented in the positive z direction, the 
vector potential is A = {B/2){—yx + xy). Expanding 
Eq. ([TJ with this choice of A yields 


H = -- 


2me 


e^B^ 
8m e 


ix^ + y^) + 


eB 

2me 


( 2 ) 


where is the 2 -component of the angular momentum 
operator, = —ih{xdy — ydx)- 

The rest of this paper uses magnetic units where length 
is expressed in units of Aq, 


Aq — 



(3) 


and energy is expressed in units of fiiOc, where oJc is the 
cyclotron frequency, Wc = eBjme. In these units, ex¬ 
pressing H in polar coordinates yields 


H = - 


1 

2 


— drXdr 
r 





(4) 


and the single-particle energy is 

i (2n -I- m -I- |m| -I-1), (5) 

where n is a nodal quantum number and m is the ro¬ 
tational quantum number about the 2 -axis. Section IHII 
shows how the Hamiltonian is modified when including 
more degrees of freedom, while Sec. II V Cl makes addi¬ 
tional modifications when expressed in hyperspherical co¬ 
ordinates. Many aspects of Eqs. and carry over 
to this formalism. 


III. N-BODY RELATIVE HAMILTONIAN 

The V-body noninteracting Hamiltonian Hpf is sepa¬ 
rable into a center of mass Hcm and relative iLrei com¬ 
ponents, 

Hn = Hcm + ( 6 ) 

where, in Cartesian coordinates akin to Eq. 



( 7 ) 
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FIG. 1. (Color online) Diagram describing the mass-scaled 
Jacobi coordinates of Eq. (H). 


Here, TVrei = TV — 1 is the number of relative Jacobi vec¬ 
tors Pj with Cartesian component s x-i and yj , and ^ is a 
dimensionless mass scaling factor , 

/ 1 \ l/.^rel 

" = («) ^ 


The center-of-mass Hamiltonian is similar in form to 
Eq. d?]), except p is replaced by TV and there is only the 
center-of-mass vector pcM- 

The linear transformation from single-particle to 
center-of-mass and relative Jacobi vectors is arbitrary, 
but in this work, the scheme used first joins identical par¬ 
ticles into pairs, then joins the center of mass of each pair 
into ever larger clusters. For odd TV, the unpaired elec¬ 
tron is joined to the center of mass of the other paired 
particles. The Jacobi vectors are labeled in a reverse 
manner, so that the last Jacobi vector, the (TV — 1)®*, is 
always the relative coordinate for a pair of particles, and 
where the Jacobi vectors of the largest clusters have the 
smallest index. For example, for five electrons in matrix 
notation this is 
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Figure [T] is a diagrammatic representation of the Jacobi 
vectors described by Eq. ®. The numbers denote parti¬ 
cle locations while the arrows denote the Jacobi vectors, 
also labeled by Pj. This choice of Jacobi tree reduces 
the size of the unsymmetrized basis needed to achieve 
antisymmetric states, as is described in Sec. IIV Bl 


IV. HYPERSPHERICAL FORM 

This section describes the hyperspherical transforma¬ 
tion of the noninteracting relative Hamiltonian, how the 
relative Hamiltonian is expressed in these coordinates, 
and the resulting adiabatic potentials. 


FIG. 2. (Color online) Semi-canonical Jacobi tree diagram. 
The sub-tree in the lower right describes the Arei — 1 nodes 
at the top of the main tree. Reading from the main node 
at R,K, any time a node is passed to the left (right), the 
coordinate picks up a factor of sinoj (cosaj). 


A. Hyperspherical coordinate transformation 

Hyperspherical coordinates are the generalization of 
spherical coordinates beyond three degrees of freedom. 
The size of the system is correlated with a single length, 
the hyperradius R, while the geometry of the system is 
encoded in the remaining degrees of freedom as a set of 
hyperangles, denoted by ft. This length R, 

Nrel 

r^ = J2pI ( 10 ) 

i=i 

is a scalar quantity and its square is the sum of the 
squared lengths of the Jacobi vectors. 

The orthogonal coordinate transformation from Carte¬ 
sian to hyperspherical coordinates has some arbitrari¬ 
ness, and many different schemes can be found in the 
literature [iTl - flH . . In this work, the semi-canonical 
construction is most useful [see also Fig. A4 of Ref. 1^1 . 
It resembles the canonical tree structure of Refs. mm, 
but instead of each new branch adding a single degree 
of freedom, each additional branch adds an additional 
two-dimensional sub-tree. In this way, the particle-like 
nature of the two-dimensional Jacobi vectors is main¬ 
tained. This set of hyperangles consists of the TV^ei az¬ 
imuthal angles (j)j associated to each Jacobi vector pj and 
the TVrei — I constructed hyperangles Oj, where 


\JeLipI 

tana, = - - 

Pj+i 


( 11 ) 


Figure [2] gives a diagrammatic representation of the 
semi-canonical construction. This Jacobi tree connects 
branches (segments) into nodes (dots), where to ev¬ 
ery node is associated a sub-length, a sub-hyperangular 
quantum number [see e.g. Eq. ®], and an angle. The 

sub-lengths Rj are defined similarly to Eq. (nnD, e.g. 
r?2 _ _i_ ^2 

R 2 — Pl+ P2- 

The Jacobi tree contains all of the information describ¬ 
ing the coordinate transformation from Cartesian to hy¬ 
perspherical. Read from the main node at {i?, TV} to 
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the top nodes at {pj, |TOj|}, any time a node is passed 
to the left (right), the Cartesian coordinate picks up 
a factor of sina^ (cosaj). The sub-tree in the lower 
right of Fig. [2] describes the N^ei — 1 nodes at the top 
of the main tree. For example, xi = piCOS(/)i. How¬ 
ever, for TV = 5, a; 1 in terms of every node is expressed 
as xi = i?sina 3 sinQ; 2 sinQ;i cos())i. The 4>j angles range 
from 0 < (j) < 2 tt while the hyperangles aj range from 
0 < aj < 7r/2. The volume element for each ipj is d(j)j, 
while for each aj it is sin^-^”^ aj cos ajdaj. 


B. Relative Hamiltonian 

The noninteracting relative Hamiltonian, Eq. o, 
transforms to 

= + + ( 12 ) 

where is the total relative z-component of the 

angular momentum, and the Laplacian operator in hy- 
perspherical coordinates ^ becomes 

= ,13) 

k is called the grand angular momentum operator [l8l| . 
Note the similarity with Eq. o, where there are radial 
and angular components. 

The eigenstates of k^ are the hyperspherical harmon- 
ics, $^^(r2), where 


and the set of hyperspherical harmonics are orthonormal 
over the hyperangles H, 

J dn $^^r(n)$^)(f2) = Sk'kSu'u. ( 15 ) 

The grand angular momentum quantum number K {K = 
0 ,1, 2 ,...) is analogous to the angular momentum quan¬ 
tum number; for the noninteracting system, it remains 
a good quantum number. Here, u is an index used 
to label the different unsymmetrized states within a 
given iV-manifold (fixed K subspace). Eor 2TVrei de¬ 
grees of freedom, there are (2TVi.ei — 2 + 2iV)(2TVi.ei -3-1- 
iV)!/(iV!(2TVi.ei—2)!) linearly independent unsymmetrized 
functions. 

The projection quantum numbers rrij associated with 
the Jacobi vectors are not good quantum numbers in the 
presence of interactions and antisymmetrization, however 
the total projection quantum number M remains a good 
quantum number of the system, even with Coulomb in¬ 
teractions (or any other interactions that depend only 
on the inter-particle distances). The set of hyperspheri¬ 
cal harmonics are made simultaneous eigenstates of 
and L™*’*°* by enforcing the constraint 

JV„i 

( 16 ) 

T=i 

\M\ < K, which is assumed in the following discussion. 
Because the center of mass has already been separated, 
it does not contribute to M, and its and energy can 
be incorporated into the system trivially. 

In the semi-canonical coupling scheme of this work (see 
Fig- HI I the unsymmetrized hyperspherical harmonics are 
expressed as 


$ 


(M) 

Ku 


(n)=n 




i=i 




n A4 




sin^'" Qffc cosl"*''+^l Ofc ^^’l™''+^l(cos2afe). 


(17) 


_ (2TVfc_|_i -|- 2k) r(nfe -|- Kk -I- |TOfe+i| -|- k) T(nk + 1) 
^ F(nfc -|- Kk + k) F(nfc -I- |TOfc_|_i| -|- 1) 


where the P are Jacobi polynomials, A4 is the normal¬ 
ization for each Jacobi polynomial, and the Kk are sub- 
hyperangular “quantum numbers,” defined recursively as 

Ki = |mi|, and (19) 

Kk = 2nk-i + Kk-i + \mk\- ( 20 ) 

The Uk [fik = 0,1,2,...) are determined after fixing the 
various Kk and \mk\. In practice, it is easier to first 
choose the Uk and |mfc|, since Uk is the order of the Jacobi 
polynomial, then determine the Kk- 


With the exception of the grand angular momentum 
K and the total azimuthal quantum number M, neither 
the Kk nor the remain good quantum numbers after 
antisymmetrizing the hyperspherical harmonics. We an¬ 
tisymmetrize the set of functions Eq. (ITT)) following the 
method of Ref. [s^ . In a closed subspace of the Hamil¬ 
tonian (here, the subspaces with both fixed K and fixed 
M), the antisymmetrized states must be linear combi¬ 
nations of the unsymmetrized basis functions within the 
same subspace. As before, the subscript u is an index 
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that distinguishes different orthogonal basis functions in 
the same manifold; each index u is associated with a 
different set of good hyperspherical quantum numbers 
that satisfy Eqs. m- Similarly, the new subscript a will 
later be used to index different antisymmetric states in 
the same manifold, but in this case, the Kk and ruk no 
longer constitute good quantum numbers. The complete 
set of basis function labels in the unsymmetrized basis 
will be indicated by a bold u = {1, 2 ,...}, and the num¬ 
ber of basis functions will be given by |m|. 

The antisymmetric basis functions can be constructed 
by first building the full matrix of the antisymmetriza- 
tion operator A = I - P 12 - P 13 - P 23 + P 123 + ■■■ 
(all A^! terms) connecting all unsymmetrized states 
in a given K,\M\ manifold. This Hermitian matrix 

Aij = is then diagonalized, i.e. we 

find the eigenvalues and their corresponding eigenvec¬ 
tors, Xia- The Afa eigenvectors with eigenvalues equal 
to N\ give the totally antisymmetric hyperangular wave 
functions, 

= E a = 1, 2, (21) 

U 

where J\fa is the number of antisymmetric states and is 
typically smaller than the total dimension of the degen¬ 
erate unsymmetrized subspace. Note that the matrix X 
depends of course on K and M, but this has been sup¬ 
pressed for notational brevity. 

The most time-consuming part of this calculation is 
the determination of the matrix Aij, which can be ac¬ 
complished using a technique proposed by Efros [s^ . In¬ 
stead of being calculated directly, the antisyinmetrization 
matrix can be found by treating Aij as many unknowns 
in a linear system of equations. The antisymmetrization 
matrix is first re-expressed in terms of a matrix equation, 

i$W(n) = ^cl,W(f 2 )^,^, ( 22 ) 

i 

where i,j € u and <i>^^(r 2 ) represents a |M|-length row 
array of all of the unsymmetrized basis functions. The 
matrix Aij is the unknown |m| x |it| matrix representation 
of the antisymmetrization operator, and Eq. (I21D evalu¬ 
ated at a single set of A^-particle coordinates constitutes 
a system of |u| equations with unknowns. The an¬ 
tisymmetrization operator is understood to act on the 
particle coordinates, so that (Jl) = (^Jl) 

Because Eq. (1221) holds true for any set of N-particle 
coordinates, substituting any random set of angular co¬ 
ordinates produces a different, linearly independent 
equation in the unknowns of the antisymmetrization ma¬ 
trix, Aij. Substituting |u| different sets of A^-particle 
coordinates into Eq. (1^^ results in a linear system of 
equations that can be solved for Aij . If we write ^jj to 
mean the two-dimensional array of unsymmetrized basis 
functions evaluated at different sets of A^-particle coordi¬ 
nates, where the column index j indexes the |m| different 


unsymmetrized basis functions, the row index 7 indexes 
the |m| different sets of coordinates, n.y, and the K and 
M quantum numbers have been suppressed, then Aij can 
be found by solving the following equation; 

( 23 ) 

i 

In practice, constructing the matrices and A^..fj is 
trivially parallelizable, but the memory requirements are 
significant since it requires the storage of two |m| x |m| 
double complex dense matrices. The size of the unsym¬ 
metrized basis increases dramatically with the number of 
particles. Choosing paired Jacobi coordinates helps limit 
the growth in spin-polarized fermion systems, since states 
with pair coordinates associated with even rrij can be 
eliminated due to symmetry. However, the growth is still 
rapid, and as a result we have only performed calculations 
in systems with up to 6 spin-polarized electrons at 1/3 
filling. In the future, with programs like SCALAPACK, 
it is feasible to push this analysis further. The problem of 
antisymmetrizing the hyperspherical functions becomes 
very challenging as unsymmetrized basis expands, and 
other strate gies for antisymmetrization may also prove 
effective IdlMdSlI. 


C. Eigenstates of the noninteracting Hre\ 

Each exact eigenfunction 'I'(i?, fJ) of the relative 
Hamiltonian, Eq. m, is separable into a hyperradial 
function, F^k{R), and one antisymmetrized hyper¬ 
spherical harmonic ,$^^(17), 

^{R,n) = (24) 

The many-dimensional hyperradial Schrodinger equation 
thus reduces to a one-dimensional uncoupled ordinary 
differential equation: 

where the noninteracting potentials (R) are given 
by 

U^k^\r)= (26) 

(X + jV,el-l/2)(JC + jV.el-3/2) ^ 2 , 1 ^ 

2Aii?2 “^8 2 ■ 

The noninteracting hyperradial solutions F^k{R) to 
the scaled Schrodinger equation, Eq. (HSl) , are 

( 27 ) 
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FIG. 3. (Color online) (a) Light solid lines show the adiabatic 
potential curves of the noninteracting four-body system whose 
minima lie below Ahwc- The lowest group {K = \M\) repre¬ 
sents the lowest Landau level. The next and next-next higher 
groups are for K = |M|-|-2 and K = |M|-|-4, respectively. The 
dark solid lines are for M = —6, the lowest of which supports 
the integer quantum Hall state, (b) The M = —6 hyperradial 
noninteracting potential curve in the lowest Landau level and 
the corresponding hyperradial energy spectrum. 

where ur (jir = 0,1,2,...) is the hyperradial quantum 
number and L is the associate Laguerre polynomial. The 
normalization, where 


Landau levels. However, these potentials are not single¬ 
particle potentials, but rather support many-body states. 
In the case of Figure |3K, energy arguments indicate the 
first excited “band” represents a four-body state with a 
filled lowest Landau level and an excited electron. 

For any K, M, the particles are confined due to the 
diamagnetic term and each of these potentials supports 
an infinite number of bound states. The noninteracting 
hyperradial solutions, Eq. (Ell), are centered within the 
non-inteacting potential wells, U^\r). For a given po¬ 
tential, each nodal excitation in the hyperradius adds one 
unit of cyclotron energy huJc, as shown in Fig. [Sjo, and 
the additional hyperradial nodes increase the overall size 
of the hyperradial wave function due to the additional 
hyperradial nodes. For a fixed M, the allowed values of 
K are K = \M\, \M\ -|- 2,... in step sizes of 2. For each 
increase in A", the confining potential moves out in the 
hyperradius, indicating the size of the system is increas¬ 
ing as the energy increases. 

Another striking feature is the potential that corre¬ 
sponds to the integer quantum Hall state. The dark 
curve with K = |M|, whose minimum is about R fa SAq, 
is isolated from the other curves with higher AT = |M| 
values. In fact, there is always a single isolated curve sep¬ 
arated from the other curves at AT = |M| = N{N — I)/2 
for all system sizes. The one antisymmetric hyperspher- 
ical wave function of this manifold, when re-expressed 
in terms of independent particle coordinates, is exactly 
identical to the wave function of the integer quantum Hall 
state. Moreover, the first excited “band”, representing a 
filled lowest Landau level and an excited electron, allows 
potentials at smaller K values. Exciting a single electron 
to a higher Landau level, allows the system to compress 
to a smaller hyperradius at the cost of one hujc unit of 
energy. 


J\f = 


urI 


r{nR + K + Nrei) ■ 


(28) 


The noninteracting many-body energies are 

= ^{2nR + M + K + N,,i). (29) 

Note that this equation is similar in form to Eq. ([5]) where 
Ur is a nodal quantum number and K plays the role of 
the \M\ term. In the limit that N,.g\ = 1, this equation 
reduces to Eq. ©. 

Figure |3] shows the effective potentials U^\r), 
Eq. (1261) . for the noninteracting system with four identi¬ 
cal fermions. The light solid lines show all the potentials 
visible in the given scale, while the dark solid lines high¬ 
light those curves with M = —6 and, from bottom to 
top, K = 6,8, and 10, respectively. One striking feature 
is that the potentials group into “bands” whose poten¬ 
tial minima group around the same energy. Each band 
is separated by a cyclotron unit of energy, reminiscent of 


D. Hyperspherical filling factor 


The most important parameter in the quantum Hall 
problem is the filling factor, r, the number of occupied 
Landau levels in a sample in the noninteracting limit. It 
is given by 


ph _ N(j)Q 
eB~ BA' 


(30) 


where p is the two-dimensional electron density, B is the 
magnetic field, N is the number of electrons, A is the 
sample area, and (jjo is the fundamental flux quantum, 
(f>o = hje in S.I. units The filling factor gives the ratio 
of the number of electrons in the system to the number 
of available single particle orbitals in the lowest Landau 
level for a given sample area. The electron density can 
be somewhat controlled with doping and gate voltages, 
but for typical experiments in gallium arsenide [l|, is 
on the order of 1 — 3 x 10^^cm“^. As an example, in a 
system with p = 2.4 x [1^, the v = 1 quantum 

Hall state is found at a magnetic field near H = lOT and 
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the V = XjZ state occurs around the much higher field 
B K, 29T. 

Experimentally, the Hall resistance quantizes to val¬ 
ues of Rh = hjve^. Defining the filling factor in terms 
of the electron density is ideal for experimental systems, 
where the local electron density is averaged over enor¬ 
mous numbers of electrons. However, for systems with 
few electrons, there is some reasonable ambiguity in es¬ 
tablishing the average density of a sample that is not 
sharply confined. In most numerical models of the pla¬ 
nar system, the ambiguity about the area is resolved by 
cutting out of the Hilbert space all single-particle wave- 
functions whose maxima in the lowest Landau level lies 
outside of a certain radius. The area of the disk defined 
by this radius is an approximate area for the few-particle 
model system. 

In the hyperspherical construction, the area of a dis¬ 
tribution of N identical mass particles correlates with 
the particles’ hyperradius, Eq. which can be used 

to define the filling factor without reference to the sin¬ 
gle particle wave functions. The connection between the 
hyperradius and the area can be established by a statisti¬ 
cal average over a large ensemble of systems, holding the 
particle number N, and the sample radius, constant. 
Each individual system in the ensemble consists of a ran¬ 
dom distribution of N particles over a disk with radius re- 
The particle density over this area is clearly p = N/nr'^, 
although the density is obviously non-uniform. 

In general, systems with different random distributions 
of N particles will have different hyperradii, but the value 
of the hyperradius squared tends to increase on average 
with increasing particle density and tends to decrease 
with disk size. The statistical average of the hyperra¬ 
dius squared for a fixed N and Trr^ over many trials is 
empirically found to be 


{R^)N,r^ 


(N - l)rg 

2 p 


(31) 


where p is the reduced mass of N particles, Eq. ([5]). For 
any individual distribution of N particles on a circle of 
radius rc, Eq. m typically does not reproduce rc accu¬ 
rately. However, the definition of the particle density is 
also only described as an average. The particle density 
p = N/ nVe for a few particles on a disk also does not ac¬ 
tually give the average density of those particles except 
with reference to a statistical average. The hyperradius 
gives a consistent way to measure an approximate area of 
a distribution of particles, even in the less clearly defined 
case of very few particles. 

Using Eq. eiD, it is straightforward to find a hyper¬ 
spherical expression for the filling factor. We first scale 
all lengths in the problem by the magnetic length of the 
system, Aq [see Eq. and get the filling factor in terms 
of the radius of a disk, that gives an approximate area 
of the sample. Substituting for this characteristic disk 


radius, equation Eq. o becomes 
N{N-1) 


(32) 


This expression for the filling factor can be used to 
distinguish the filling factors of different noninteracting 
wave functions. In the lowest Landau level, the (unnor¬ 
malized) noninteracting hyperradial wave functions have 
no hyperradial excitations and depend very simply on the 
hyperradius and the grand angular momentum: 

+ oc (33) 


The peak of this wave function is easily located at = 
2K/p, and indicates a most-likely hyperradius for each 
noninteracting wave function. Substituting this R^ into 
Eq. (15^ gives the following simple expression for the 
filling factor in the lowest Landau level. 


N{N-1) 

2if 


(34) 


This formula accurately determines the filling factor for 
the 1 / = 1 and the Laughlin z/ = 1/m states, for both bo¬ 
son and fermion systems. Other Jain composite fermion 
fillings are not accurately assigned by this function; in¬ 
stead the identified states occur at a small shift away 
from this ideal hyperspherical filling function (see Ap¬ 
pendix 1^. This shift is most prominent in few-body sys¬ 
tems, and approaches zero in the thermodynamic limit. 


E. Coulomb interaction 


Including Coulomb interactions in the noninteracting 
Hamiltonian Eq. (IT^ yields 

where k is a dimensionless parameter that determines the 
strength of the Coulomb interactions. Here, 


K — ^ — °o Rh 1 _ 1 j Eh 

AneXo fiuic XofujJc^/^o e/coV ^c’ 

where oq is the Bohr radius, Eh is the Hartree unit of 
energy, and o^Eh = XqHuJc = ^^/(2me) with QoEh = 
e^/(47reo). e is the permittivity of the material, in units 
of the permittivity of free space eg. In the lowest Landau 
level for typical experiments in gallium arsenide, n is on 
the order of one or smaller; for example, k rs 0.76 at 
B = lOT. 

The form of the hyperspherical Coulomb term (7(0) 
depends on the choice of Jacobi vectors and hyperangles. 
In general, from single-particle to relative coordinates the 
Coulomb interaction involves a linear combination of only 
the relative Jacobi vectors. 


E 





E 


K 


(37) 












For a concrete example, invert the transformation matrix 
Eq. ([H]) to express the j = 1... 5 in terms of the Pk, 
k = 1.. .4: and take vector differences. 

Because we use the antisymmetrized hyperangular 
functions the Coulomb interaction need only be calcu¬ 
lated between one pair of electrons, then scaled by the 
number of pairs. It is simplest to use the last Jacobi vec¬ 
tor that is proportional to the distance between a pair of 
electrons. In the reverse construction of the hyperangles, 
the last Jacobi vector is the simplest to express in hy- 
perspherical coordinates [see Eq. m and Fig. [2], such 
that 

^ ^_, (38) 

R 2 V 2/r i?cosaAr„j_i 

This simplification is valid only when the basis fnnctions 
are either totally symmetric or totally antisymmetric. 
Thus, integrating the above expression in the basis of 
unsymmetrized hyperspherical harmonics [see Eq. (13] 
reduces to a one-dimensional integral in because 

the other integrations are accomplished via orthogonal¬ 
ity of the angular functions. In practice, Gauss-Jacobi 
quadrature is used to evaluate the integral in daN„i-i- 

The strategy to diagonalize Eq. (IMll remains the same 
as that for the noninteracting system. First, M remains a 
good quantum number and each M block of i/pei is diag¬ 
onalized independently. However, the expansion Eq. (1241) 
is no longer strictly separable into radial and hyperangu¬ 
lar functions. Instead, the hyperangular channel func¬ 
tions depend parametrically on R, where 

4f{R, n) = i?-^r,l-Sl /2 (gg) 

X 

Here, y labels each channel and the channel fnnctions 
(i?; 17) are orthonormal for a fixed hyperradius, 

J dn 17)$^f) {R; 17) = . (40) 


The hyperradius R is treated as an adiabatic parame¬ 
ter, where the adiabatic Hamiltonian iJad, 


i?ad = 


2pR‘^ 


K^ + (N,,i-l/2)iN,ei-3/2)\ (41) 


d- r,2 I 1 A ,r I ^(17) 


+ ^R^ + -M + K 


R 


is diagonalized at each value of the hyperradius. To find 
eigenstates of Eq. gD), the channel fnnctions are ex¬ 
panded at a fixed hyperradius using the antisymmetrized 
hyperspherical harmonics, 


$W(i?;17) =Y,CKaW^KJm- (42) 

Ka 


Under this expansion, the matrix elements of i^ad 


= + (43) 

H 




R/Ao 


FIG. 4. (Color online) Top: Adiabatic potentials Ux (R) 
for N = S, M = —9 and k — 1. Bottom: Adiabatic potentials 
Ux~^^\R) for N = 4, M = —18 and k = 1. Similar to 
Fig- El the larger gaps (on the order of huuc) indicate magnetic 
excitations, while the smaller splittings are due to Coulomb 
interactions. The separate clusters of curves (distinguished by 
color online) indicate different K manifolds, with K = \M\ for 
the lowest grouping of curves and K increasing from bottom 
to top. 


The brackets indicate that the integrals are taken only 
over the hyperangles and the a’s that enumerate the dif¬ 
ferent hyperspherical harmonics of a given K manifold 
are distinguished by a subscript. 

The solutions at each R describe a set of adiabatic 
potentials U^\r), where 


= Ul^^\R)^‘f'>{R; n). (44) 


Fig. HKa) shows the adiabatic potentials U^\r) for 
N = 3, M = —9, and k = 1, while Fig SKb) shows the 
adiabatic potentials for iV = 4, M = —18, and k = 1. 
The interactions are weak such that the different K man¬ 
ifolds are still distinguishable. However, the Coulomb 
interaction has split the degeneracy of the potentials. If 
K were increased, then the states comprising different K 
manifolds would begin to overlap. 

Without approximations, the many-dimensional 
Schrodinger equation with Coulomb interactions reduces 
to an infinite set of one-dimensional coupled ordinary 
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differential equations in terms of the adiabatic potentials: 


- ^ E {2 ^xx' ^ + Qxx' } 4 ? (^) = 0 . (45) 

^ x' ^ J 


The T’xx' Qxx' niatrices, 


Pxx'(R) 



^ \ 

ai? / 


(46) 


and 


*3xx' ('^) 



ai?2 /’ 


(47) 


describe the non-adiabatic coupling between different 
channel functions, the 

A good approximation to the adiabatic potentials is to 
neglect the coupling between different K manifolds and 
apply degenerate perturbation theory. Diagonalizing the 
Coulomb matrix in each manifold, with restricted matrix 
elements {Ka'\C{n)\Ka), yields the eigenvalues C^\ 
where the subscript 7 labels the different eigenstates of 
the Coulomb matrix in a fixed K manifold. The Coulomb 
interaction under this approximation reduces the prob¬ 
lem to a set of uncoupled one-dimensional potentials, just 
like in the noninteracting system. The resulting adiabatic 
potentials are 


/ /^{^) \ 

(i?) Skk' f (R) + j . (48) 

Fig. [S{a) shows the lowest hyperangular Coulomb 
eigenvalues, the , for TV = 6 particles obtained by 

diagonalizing the Coulomb interaction within the degen¬ 
erate manifolds having K = —M. The minimum eigen¬ 
values are shown from the iV = 15 manifold (filling factor 
iz = 1) to the TV = 45 (jz = 1/3). The values shown 
are dimensionless quantities. A classical minimization of 
the Coulomb potential at fixed hyperradius gives a lower 
bound to the quantum eigenvalues, and the following for¬ 
mula gives approximate minimum values, namely: 

C„ii„(TV) « (0.12 -k 0.33/Af)Af2(TV - 1) (49) 

For instance, for TV = 6 ( 8 ) particles, the direct minimiza¬ 
tion gives Cmin = 31.7824(71.5427), whereas the approx¬ 
imate formula gives 31.5 (72.2). The value of 
computed with hyperangular wavefunctions is within 7% 
of this minimum value. The classical minimum hyper¬ 
angular Coulomb potential corresponds to the internal 
geometrical distribution of particles that minimizes the 
energy for any size system. As \M\ is increased, the in¬ 
creasing Hilbert space size allows the particles to better 
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FIG. 5. (Color online) (a) The minimum hyperangular 
Coulomb eigenvalues, the dimensionless in Eq. (148 II for 

N = 6 particles, obtained by diagonalizing the Coulomb in¬ 
teraction within the degenerate manifolds having K — —M 
(which corresponds to the lowest Landau level). Coulomb 
eigenvalues are shown for the K = 15,..., 45 manifolds, which 
correspond to filling factors v = 1... 1/3. (b) Two Yrast 

plots, plots of the minimum energy at each M, of the A, M 
manifolds for six particles in the lowest Landau level using 
different approximations. The energies of the upper (black) 
curve are calculated entirely in first order perturbation theory 
in the lowest Landau level. For the lower energies (red), the 
hyperangular energies are calculated perturbatively (taken 
from (a) above), but the hyperradial energies are calculated 
using exact numerical techniques. The lower yrast spectrum 
constitutes a departure from the strict lowest Landau level 
approximation. In this model, no confinement potential is in¬ 
cluded, and particle confinement occurs only due to angular 
momentum conservation. 


approach this ideal, energy minimizing internal geome¬ 
try. 

Including the potentials of Eq. (|4^ in Fig. |4] only 
slightly alters the potentials, and the changes are largest 
for the higher energy potential curves within a given K 
manifold. The hyperangular eigenvector corresponding 
to the smallest eigenvalue is the state that minimizes 
the Coulomb interactions. For example, we find 98% 
overlap with the hyperangular numerical ground state 
and the hyperangular part of the Laughlin function for 
the three-body system with K = —M = 9, indicating 
that the ground state in this system is a quantum liquid. 
The energies the A, M manifold restricted calculation 
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can be found by multiplying the hyperangular eigenval¬ 
ues, , times the hyperradials energies, which can 

be found in different ways. Ignoring the changes to the 
potential curves due to the Coulomb interactions consti¬ 
tutes treating the hyperradial energies to first order in 
perturbation theory. 

In the lowest Landau level, hyperspherical energy cal¬ 
culations using first order perturbation theory match 
the results of conventional planar configuration inter¬ 
action calculations, which are considered the numerical 
standard in quantum Hall studies. In the hyperspheri¬ 
cal treatment, perturbation theory calculations consists 
of multiplying the hyperangular Coulomb eigenvalues, 
times the first order perturbation calculations of 
the hyperradial expectation value of 1/r from Eq. 1481 
The upper (black) Yrast plot in Fig. HKb) for a six parti¬ 
cle hyperspherical system under these approximations is 
identical to the exact numerical diagonalization plot in 
Fig.3 of Ref. to within the numerical accuracy shown 
in their Table III [1^ . We present a few numerical exam¬ 
ples for selected systems in the third row of TableH) These 
values are identical to within published numerical accu¬ 
racy to standard planar numerical calculations [U 13 > 
given on row 2, and are significantly more accurate than 
calculations performed using a spherical geometry, given 
in row 1. Numerical studies in the spherical geometry of 
Ref. @ are another standard in the quantum Hall sys¬ 
tem, but while calculations on the Haldane sphere cap¬ 
ture much of the physics of the quantum Hall system, 
they are not considered numerically accurate due to fi¬ 
nite size effects. The energy shifts given in row 1 are 
extrapolations from the spherical geometry to an infinite 
plane, performed using leas t sq uares quartic fits to the 
energy shifts following Ref. [3l- Despite the numerical 
inaccuracy of the spherical geometry, we present these en¬ 
ergy shifts as another comparison to standard techniques 
in the quantum Hall field. 

The calculations can be improved from the lowest Lan¬ 
dau level restriction by treating the hyperradial ener¬ 
gies using exact numerical techniques. In this approx¬ 
imation, we maintain the single K, M manifold degener¬ 
ate perturbation theory approximation, but consider how 
the resulting Coulomb hyperangular eigenvalues alter the 
hyperradial potential curves. Solving the hyperradial 
Schrodinger equation using exact numerical techniques 
yields a slightly lower energy than the energies given us¬ 
ing pure perturbation theory, as shown in both Fig. IlKb) 
and row 4 of Table ID This technique constitutes incor¬ 
porating additional functional space from higher Landau 
levels into the Hilbert space of the problem. This ap¬ 
proximation would be analogous to a variational energy 
minimization based on a scaling parameter multiplying 
all of the single particle radial coordinates from the origi¬ 
nal Slater determinant basis. The energies are lowered in 
this approximation because the introduction of Coulomb 
repulsion expands the hyperradial potential curves out¬ 
ward, and the resulting hyperradial energies are lowered 
by this expansion. The solutions are still confined by the 


diamagnetic term of the Hamiltonian, but no additional 
confinement potential has been included to compress the 
electrons into a smaller area. Including a confinement 
potential to model the inward Coulomb pressure due to 
the bulk might be appropriate in modeling many body 
condensed matter systems, but it is less appropriate for 
calculations in trapped atom systems or quantum dots. 

More accurate calculations require extending beyond 
the single K, M manifold approximation, and are not 
trivial. However, upper and lower bounds on the exact 
energies of the system can be established using additional 
approximations. A lower bound can be established using 
the Born-Oppenheimer approximation [3 - In this ap¬ 
proximation, higher K manifolds with the same M are 
included in the hyperangular energy calculations, and the 
hyperradial energies are calculated by solving Eq. US] nu¬ 
merically while neglecting the nonadiabatic coupling P 
and Q matrices. If the fixed-i? hyperangular calcula¬ 
tion is fully converged, these values can be proven to 
be lower bounds on the ground state energies for each 
of the corresponding symmetries. A few example lower 
bounds are given in row 5 of Table |T| The upper bounds 
given shown in row 6 of the same table are calculated 
under the adiabatic approximation, which differs from 
the Born-Oppenheimer approximation by including the 
diagonal elements of the Q matrix in the numerical solu¬ 
tion to Eq. US] When fully converged, these values give 
strict upper bounds to the ground state energy for those 
symmetries. As before, these upper and lower bounds 
apply to a system with no additional confining potential 
besides the diamagnetic confinement. 

The fact that the upper and lower bounds in this hy¬ 
perspherical representation differ by less than 10“^ for 
all cases shown is strong evidence that the adiabatic rep¬ 
resentation is unusually effective in this system. By com¬ 
parison, the ground state energies of the hydrogen nega¬ 
tive ion were found to differ by 1.7% in the corresponding 
upper and lower bound levels of the approximation, (see 
j51|i. A major reason why the adiabatic formulation is 
much more accurate for the quantum Hall problem is that 
the charged particles here do not experience any attrac¬ 
tive potentials, which cause potential valleys in the H“ 
system that have much stronger nonadiabaticity. 

Another indication that the adiabatic hyperspherical 
representation is particularly appropriate is that the adi¬ 
abatic potential curves show only weak nonadiabatic cou¬ 
pling between different manifolds, as evidenced by a rel¬ 
ative lack of avoided crossings. The weak coupling be¬ 
tween different manifolds makes it difficult to discern any 
avoided crossings in Fig. 0] on the scale shown. To better 
visualize any avoided crossings. Fig. [6] shows the scaled 
adiabatic potentials g^\R), 

(R) = R^ ([/W (i?) - [{g/8)R^ + M/2]) , (50) 

for the same three- and four-body systems as in Fig. 0] 
(M = —9 and M = —18, respectively, each with k = 1), 
though on a much larger scale in R. At R = 0, the 
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TABLE I. Comparisons to conventional systems. All rows are estimations of the ground state energy change due to the 
introduction of the Coulomb interaction,converted to our system of magnetic units, defined in the text. Row 1:Extrapolations 
of AE from Haldane sphere configuration interaction calculations, using quartic, least squares fits following the method in 
Ref. [^. The extrapolations become less accurate as the computational size of the systems limit the number of different 
sphere sizes that can be calculated. Row 2: Planar exact numerical diagonalization AE in the lowest Landau level: first 
two values from Ref. [47|*f3/\/2i. and last three values from Ref. [i^. Row 3: Degenerate perturbation theory calculation 
using hyperspherical picture within a fixed-K manifold along with a first order perturbation treatment of the hyperradial 
equation. Row 4: Degenerate perturbation theory within the fixed-K manifold with an exact, nonperturbative treatment of the 
hyperradial differential equation. Row 5: Born-Oppenheimer approximation neglecting nonadiabatic coupling matrices P and Q 
with exact non-perturbative treatment of the hyperradial differential equation, which constitutes a lower bound approximation 
when well-converged. Row 6: Adiabatic approximation, which includes the diagonal element of the Q matrix in the lowest 
potential curve, with an exact, nonperturbative solution of the hyperradial differential equation. These constitute strict upper 
bounds to the ground state energies. 


N,-M 

3,9 

3,15 

4,18 

5,30 

6,45 

AE, Haldane sphere, fit, extrapolation 

0.71656 

0.5526 

1.310 

2.04 

« 3 

AE, Planar calculations^ 

0.716527 

0.55248 

1.30573 

2.02725 

2.86015 

AE, Perturbation Theory 

0.716527 

0.55248 

1.30573 

2.02725 

2.86015 

AE, Degenerate fixed-E 

0.704637 

0.54792 

1.28552 

1.99742 

2.81994 

AE, Born-Oppenheimer (lower bound*) 

0.70198 

0.54722 

1.28086 

1.99226^ 

- 

AE, Adiabatic (upper bound) 

0.70204 

0.54723 

1.28092 

1.99230 

- 


® First two values are (3/’\/2)times values taken from |47ll . the remaining 3 values are taken from [S 
The value shown may not be a converged lower bound to all digits shown. 


system reduces to the eigenvalues of the operator 
[see Eq. ([M]) ]. while at small R the gx{R) are linear in R 
with slopes given by . All of the are positive, 

so the curves in Fig. [6] have been shifted by the smallest 
eigenvalue of the lowest K manifold to put the curves 
within the same scale. 

With this scaling, many close avoided crossings be¬ 
come visible through the higher K manifolds. In gen¬ 
eral, most of the crossings appear diabatic in nature. 
These crossings occur at much larger hyperradii (here, 
kR > 50Ao) than the scale of Fig. [4] (i? < llAo). For 
comparison, 99% of the noninteracting hyperradial wave 
function is contained within R < 8.35Ao for 3 particles 
and R < 10.25Ao for 4 particles. Even well outside of this 
region, the lowest curve of the lowest K manifold remains 
isolated from the rest, which suggests the adiabatic ap¬ 
proximation is a good approximation for the ground state 
of the lowest Landau level. In addition, although Fig. [6] 
has K = 1, the curves are universal, and their behavior 
changes smoothly with a trivial scaling in k. In other 
words, the excellence of the adiabatic approximation ap¬ 
plies within the relevant potential range for a very wide 
range of magnetic fields. 

A quantitative measure of the non-adiabatic 
coupling strength is the dimensionless quantity 
P^^,{R)/[ui^\R) - U^f\R)]. Figure [7] shows this 
quantity as a function of kR for N = 3 and M = — 9, 
and N = 4 and M = —18 for y = 1 and y' = 2,... , 6, 
that is, the five curves represent the coupling strength 
from the Laughlin 1/3 potential (the lowest channel) to 
the five next lowest channels. Like the potentials shown 
in Fig.ini these curves are also universal in the sense that 


they scale simply as a function of k. Even the strongest 
nonadiabatic coupling is very small compared to unity, 
and the coupling to other channels is even weaker still, 
indicating the validity of the adiabatic approximation. 
Other higher x' are not shown as their coupling strength 
is weaker than 0(10“®). 

To conclude this section, the Coulomb interaction acts 
to split the states within a given K manifold, yet does not 
lead to strong nonadiabatic coupling within the region 
where the potentials are deepest. Even with Coulomb 
interactions, the separation of hyperradial and hyperan- 
gular degrees of freedom is an excellent approximation. 


V. EXCEPTIONAL DEGENERACY 

One of the benefits in describing this system in hyper¬ 
spherical coordinates is that the set of antisymmetrized 
hyperspherical harmonic basis functions in any K, M 
manifold forms a complete basis in the absence of interac¬ 
tions. From perturbation theory, it is well known that, in 
a set of functions, turning on interactions will typically 
act to lower the energy of the ground state relative to 
all the higher energy states. This effect is strengthened 
by the presence of additional degeneracy in the system. 
If the basis functions prior to turning on interactions are 
degenerate, then increased degeneracy in the noninteract¬ 
ing picture should lead to an increased energy separation 
of the ground state. In other words, K, M manifolds with 
enhanced degeneracy relative to their neighboring K, M' 
manifolds should exhibit more strongly gapped ground 
states. As a result, we predict that manifolds with excep¬ 
tionally high degeneracy are likely to also be identifiable 
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FIG. 6. (Color online) Top: N = 3 scaled adiabatic po¬ 
tentials gx{R) [see Eq. (1501) 1 shifted by the smallest eigen¬ 
value Ci_i from diagonalizing the Coulomb interaction in the 
K = —M = 9 manifold. Bottom: = 4 scaled adiabatic 

potentials gx{R) shifted by the smallest eigenvalue Gi®! from 
diagonalizing the Coulomb interaction in the K = |M| = 18 
manifold. Alternating thick (red) and thin (blue) lines label 
different K manifolds. 


quantum Hall states. 


A. Exceptional degeneracy derivation 

The following details how the exceptional degeneracy 
is derived starting from group theory. Only the lowest 
Landau level is considered in this paper, that is, only 
those states with K = —M. 

First, we derive the discrete function of |M| that de¬ 
scribes the growth in the number of antisymmetric states. 
These integer sequences are intimately related to generat¬ 
ing functions. For example, from combinatorial consid¬ 
erations, the generating function Gn{x) for the overall 
degeneracy for a fixed N of spin polarized fermions in 
the lowest Landau level can be derived using integer par¬ 
titions. In the lowest Landau level, each unsymmetrized 
iC, M manifold of the relative hyperangular functions 
times forms a complete, translationally invariant ba¬ 
sis of polynomials in N = N^ei variables (p * e*'^) that 
are homogeneous of order K = |M|. According to Eq. 
(25) of |52l |. the degeneracy of the symmetric irreducible 
representation of this basis is equal to 

dsym{.K,N)=pN{K)-pN{,K-l), (51) 

where pN{k) is the number of partitions of the integer K 
into parts no longer than N. The number of partitions 
can be calculated using a generating function (^ . 



kRIAq 

FIG. 7. (Color online) P^^,{R)/[U^’^\R) - f/^f’(i?)] as a 
function of kR for N = 3 and M = —9 (top) and N = 4 and 
M = —18 (bottom) for y = 1. Solid, long-dashed, dotted, 
dash-dotted, and short-dashed lines are for y' = 2, 3,4, 5 and 
6, respectively. 


-| OO 

Zn{x) = x^Pn{K). (52) 

i=i ^ ^ K=i 

Combining dSU and (l5^ with the fact that there is a 
one-to-one mapping between the symmetric irreducible 
representation at K and the antisymmetric irreducible 
representation at A -|- N{N — l)/2 yields the generating 
function 

(53) 

i=2 ^ 

for spin-polarized fermions in the lowest Landau level. 
The coefficients of the Taylor sequence of the generat¬ 
ing function yields the integer sequence function whose 
elements are equal to the number of degenerate an¬ 
tisymmetric states in any given K = \M\ manifold. The 
coefficients are equivalent to 77^ defined in section 
IIVBI The way this generating function is used is to ex¬ 
pand the above product in powers of x, namely 

CO 

Gn{x)= (54) 

|M|=0 

We have verified by brute force computation that the re¬ 
sulting integer coefficients are precisely equal to the 
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FIG. 8. (Color online) Total degeneracy of antisymmetric 
states for the four-body system in the lowest Landau level as a 
function of \M\. Solid and dashed lines show upper and lower 
envelope functions, respectively, while the points show the 
number of degenerate antisymmetric states in that manifold 
for each value oi \ M\ = K. 


degeneracy for many values of N and \M\. The degenera¬ 
cies and generating functions can alternately be Yr de¬ 
rived using group theory [H^. The points of Fig. [8] show 
the number of antisymmetric states for the four-body 
system as a function of \M\. The first non-zero value at 
\M\ = 6 is the integer quantum Hall state. There is only 
the one antisymmetric hyperspherical harmonic function 
with K = —M = 6 as is expected for a closed shell in the 
independent electron picture. Two other notable points 
on this scale include those at \M\ = 18 and 30, which 
corresponds to the Laughlin 1/3 and 1/5 states, respec¬ 
tively. Otherwise, the general trend as \M\ increases is 
that the number of antisymmetric states oscillates about 
an overall polynomial growth. 

There are many approaches to quantify the small vari¬ 
ations in degeneracy on top of this polynomial growth, 
such as comparing the degeneracy of nearest-neighbors or 
making comparisons after dividing out the largest power 
in \M\. We choose to derive two polynomial functions 
that envelop the degeneracies, and then compare the rel¬ 
ative heights above the lower envelope. The top envelope 
function is forced to go through the integer quantum 
Hall point at \M\iqh = N[N — l)/2, while the bottom 
function is forced to go through the zero degener¬ 
acy value at \M\iqh + 1. The solid and dashed lines 
of Fig. IS] show the upper and lower envelope functions, 
respectively [see also Eqs. (|57)l and (158)) ]. 

The envelope functions are derived directly from 
the exact integer sequence degeneracy function 
a|^|. For small systems, Mathematica can usu¬ 
ally find the sequence functions directly by using 
SeriesCoefficient[G'Ar(a;), {x, 0, |M|}] (—M is assumed 

to be non-negative). In our experience, it is easier to 
first do a partial fraction decomposition of Eq. (1551) 


using Mathematica’s Apart [Gat (cc)], and then find the 
series coefficient of each term and add them together. 
Regardless, the sequence function aj'^j is grouped in 
powers of |M|. The coefficients of the upper (lower) 
envelope polynomial are derived by evaluating the 
coefficients of a[^| at \M\iqh {\M\iqh + 1). 

As a concrete example, for four particles the partial 
fraction decomposition of G/i(x) is 


Gi{x) 


1/24 

(1-X)3 ■ 
, -5/32 
1 -I- X 


-1/8 23/288 1/16 

(1 — xY 1 — X (1 -b xY 
-(1-3^) , 2 + x 

8(l + x2) ^ 9(l + x + x2)- ^ ^ 


Converting Eq. (l55l) to a|^| (using Mathematica func¬ 
tions like ExpToTrig) yields 

m2 (_i)|a^I_i 


j(4) _ 

|M| - 48 


36 sin 


r|M| 


16 

— 36 cos 


|A^| + ^(-27(-1)I"^I-1 

TrjMj 


64 cos 


(56) 




The upper envelope function comes from evaluating 
the coefficients of Eq. (1561) at |M|iqh. Here |M|iqh = 6, 
such that the coefficient of the |M|2 term is a constant 
1/48 and the coefficient of the \M\^ term is 0. The co¬ 
efficient of the |M|° term is determined last. It is found 
by forcing the polynomial to equal 1 when evaluated at 
|A1|iqh, specifically 36/48 -l-0-l-x = l, soa value of 
X = 1/4 is the final term. This yields an upper envelope 
function of 


Vi “ 


|M|2 

48 


1 

4' 


(57) 


The lower envelope function comes from evalu¬ 
ating the coefficients of Eq. (1551) at |M|iqh + 1. Here 
|A1|iqh + 1 = 7, such that the coefficient of the |Mp term 
is a constant 1/48 and the coefficient of the |M|^ term 
is —1/8. The coefficient of the \M\^ term is determined 
last. It is found by forcing the polynomial to equal 0 when 
evaluated at |M|iqh + 1, specifically 49/48 —7/8-1-x = 0, 
so a value of x = —7/48 is the final term. This yields a 
lower envelope function of 


(4) |Mp |M| 7 

-|M| 48 8 48' 


(58) 


Defining the upper envelope function as and 

the lower envelope function as L]^{M)^ we have derived 
the following expressions in terms of the coefficient Gn = 
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N{N — 1)/{NI)‘^ of the maximum power \M\^ ^ as: 

U 3 /C 3 = |M| + 3 

L 3 /C 3 = |M| - 4 

U 4 /C 4 = + 12 

L 4 /C 4 = - 6 |M| - 7 

U 5 /C 5 = |M|3 - 9M^ + 36|M| + 260 

L 5 /C 5 = |M|3 - 9M^ - 9\M\ - 143 

Ue/Ce = M^- 20 |Mp+ ISOM^- 180|M|+ 3105 

Lq/Cq = M^- 20|M|3 + 60M2 + 80|M| - 256 

U 7 /C 7 = \M\^ - 75M^/2 + 1340|M|V3 - 825M^ 

- 5865|M| + 328293/2 

L 7 IC 7 = |M|® - 7hM‘^l2 + 1340|M|3/3 - 2400M2 
+ 6560|M| - 323576/3 

t/g/Cs = M® - 63|M|® + 1400M^ - 10920|M|® 

- 21168M2 + 784000|M| + 5234432 
Lg/Cg = M® - 63|M|® + 1400M^ - 14070|M|® 

+ 78057M2 - 552475|M| + 714734 
U 12 /C 12 = - 275|M|® + 63195M®/2 

- 1960200|M|^ + 70516292M® 

- 1440719280|M|® + 13562493120M^ 

+ 40317868800|M|® - 2246025672000M2 
+ 20132954569728|M| + 289846790411904 
L 12 /C 12 = - 275|M|® + 63195M®/2 

- 1960200|M|^ + 70516292M® 

- 1460365830|M|® + 16263893745M^ 

- 99340265200|M|® + 962578332375M2 

- 11917386104239|M| - 101761341423733/2 


From the envelope functions, the relative degeneracy 
g^ei can now be defined as 


ffrel = 


AN) _ JN) 
^\M\ 


^|M| 


= W JN) 


‘'\M\ 


^\M\ 


(59) 


that is, the relative height of above the lower en¬ 
velope, with respect to the separation between the two 
envelope functions. 


B. Connection between degeneracy and energy 
gaps 

Figure [5] illustrates the connection between the en¬ 
ergy gaps that appear when solving Schrodinger’s equa¬ 
tion for different M and the exceptional degeneracies. 
Panel (a) shows the energies from solving the four-body 
Schrddinger equation in distinct iF, M manifolds for k, = 
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FIG. 9. (Color online) Connection between exceptional de¬ 
generacy and energy gaps for N = 4 electrons, (a) The dashes 
indicate the energies at each \M\. The IQH, 1/3 and 1/5 fill¬ 
ing fractions are labeled, (b) The energy gap is defined as the 
smallest energy difference between the ground state at \M\ 
and all other states with \M'\ < \M\. This could be neg¬ 
ative if the comparative state is lower in energy. The IQH 
state’s energy gap is not shown, but is about IfkOc, since the 
next nearest state is approximately one vibrational quantum 
away, (c) The relative degeneracies are defined in Sec. IV Al 
see specifically Eq. (ESJ- 


1 and various \M\, but ignoring coupling between K man¬ 
ifolds. The IQH state is recognizable at |M| = 6 as it is 
non-degenerate and isolated from the rest. The Laughlin 
1/3 state at |M| = 18 also shows a significant lowering 
of the energy compared to neighboring M values. Panel 
(b) shows the energy gaps, that is, the smallest energy 
difference between the ground state at |M| and all other 
states with \M'\ < \M\. The IQH state is not shown 
on this scale as its energy gap is ~ hwc- Again, the 1/3 
state, among others, shows a prominent positive energy 
gap. Panel (c) shows the relative degeneracies. The IQH, 
1/3, and 1/5 states all have unity relative degeneracy, 
while other states also show prominently. Remarkably, 
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FIG. 10. (Color online) Relative degeneracies [see Eq. (I59II ] for 
the six-body system are shown as a function of |M|. Squares 
show the integer quantum Hall effect and the Laughlin v — 
1/3,1/5, ... states; triangles show the Jain states of two filled 
composite Landau levels (also called A levels); circles show 
the remaining unidentified states. 


although the upper envelope is derived based only on the 
IQH state, it passes through all Laughlin-type degener¬ 
acy points as well, this precise matching of the upper 
envelope function does not hold for all higher particle 
numbers N. In general, we designate those K = —M 
states with relative degeneracy close to one as the ex¬ 
ceptionally degenerate states, and suggest that these are 
candidates for observable iV-body ground states that are 
markedly lower in energy than their neighbors. 

Figure [TUI shows the relative degeneracy as calculated 
from Eq. (lUUIl for the six-body system as a function of 
|M|. The squares identify the degeneracies in systems 
with values of M that include the integer quantum 
hall and Laughlin states. The triangles identify the 
degeneracies of systems that include the Jain states 
of two filled composite fermion Landau levels, which 
also show a large relative degeneracy as compared to 
their neighbors (which appear near g^ei Ri 0) and their 
next-nearest neighbors. Circles show the remaining 
unidentified states. A brief discussion of the identifi¬ 
cation of Laughlin and Jain states is given in appendix lAl 


VI. CONCLUSIONS 

For testably small systems in the lowest Landau level, 
many of the K, M manifolds of antisymmetrized func¬ 
tions contain the identifiable Laughlin and Jain com¬ 
posite fermion states for few-body systems. The K, M 
manifolds of the identifiable Laughlin and Jain composite 
fermion states all exhibit exceptionally high degeneracy 


compared to the majority of other K, M manifolds in the 
lowest Landau levels. Although the relative degeneracy 
does not uniquely identify the known composite fermion 
filling states, our results suggest that high degeneracy 
plays a role in strengthening the energy gap of observed 
and described fractionally filled states. As such, it may 
also be interesting to examine the low-lying ground states 
that are associated with exceptionally high relative de¬ 
generacy, but which are not associated with the Laughlin 
or Jain sequences. 

The adiabatic hyperspherical potential curves we have 
calculated are astonishingly devoid of strong couplings 
and avoided crossings, which is a sign that the adiabatic 
approximation is extremely and unusually accurate for 
this system in the parameter range typically probed in 
FQH experiments. In other words, the hyperradial de¬ 
gree of freedom is accurately quasi-separable in the exper¬ 
imentally studied range of the FQHE. Although the hy¬ 
perspherical adiabatic approximation in the lowest Lan¬ 
dau level reproduces the results of exact numerical config¬ 
uration interaction calculations, the treatment accesses 
different aspects of the problem than either the Slater 
determinant construction or the composite fermion pic¬ 
ture. The quasi-separability of the hyperradial degree of 
freedom is a feature not considered in other treatments 
of the quantum Hall system, and the hyperradius does 
not naturally emerge from either the Slater determinant 
or the composite fermion constructions of the quantum 
Hall system. At this time, we cannot draw any clear 
connections between our construction and the composite 
fermion framework, but suggest that this alternate per¬ 
spective on the problem may allow us to examine prop¬ 
erties of the quantum Hall system that emerge more nat¬ 
urally out of the hyperspherical picture. In particular, 
the hyperspherical construction suggests the existence of 
higher energy states that are hyperradial excitations of 
the ground state wave functions. These hyperradially 
excited states should have the same internal structure 
as their ground-state counterparts; excitations between 
such states should represent a breathing mode that could 
be probed spectroscopically. These aspects will be ex¬ 
plored in future publications. 
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Appendix A: Relative M and Identification of 
Quantum Hall States 

The identification of the experimentally observed frac¬ 
tional quantum Hall states in systems with a modest 
number of particles is not trivial. Although the high 
exceptional degeneracy of a K,\M\ manifold is highly 
correlated with the presence of a quantum Hall ground 
state, it is not demonstrated to be a diagnostic of the 
presence of a quantum Hall state. In addition, the filling 
fraction as given by Eq. (1301) is correct in the thermody¬ 
namic limit but is only approximately correct for small 
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TABLE II. Sample list of identified N-body quantum Hall 
states in the lowest Landau level. M is the total relative 
azimuthal quantum number of Laughlin and Jain states iden¬ 
tified by exact numerical diagonalization in a spherical geom¬ 
etry [^. ucF gives the filling factor of identified QH states 
according to the Jain composite fermion picture, including a 
correction that accounts for the finite size shift associated with 
the spherical geometry, vhs is the calculated hyperspherical 
filling factor, given by Eq. The final column gives a finite 
size correction to the hyperspherical filling factor. 



TABLE III. Jain Composite Fermion states in the lowest Lan¬ 
dau level. The M for an A-particle system at filling fraction 
iz = j/*/(l + 2pv*) are given for the composite fermion states 
most strongly observed in experiments, where p is an inte¬ 
ger indicating the number of pairs of composite fermion flux 
tubes attached to each electron. 


V* 

-M 

Restrictions 

1 

^(^-i)(2p+i) 


2 

A(^+p(A-l)) 

even N only 

-2 

Ni-^+p{N-l)) 

even N only 

3 

A(^+p(A-l)) 

N mod 3 = 0 only 

-3 


N mod 3 = 0 only 


systems. It is also of limited use for uniquely identifying 
the quantum Hall ground states. Instead, the fractional 
quantum Hall states of important filling factors are iden¬ 
tified by using results from conventional, exact numerical 
diagonalizations in finite systems using planar, spherical, 
or toroidal geometry. 

For example, in a system of 6 particles, Eq. dSO]) 
would predict that the v = 1/3 state should appear 
when the single particle Hilbert space is restricted to 18 
orbitals in the lowest Landau level, or, in other words, 
when the number of magnetic flux quanta in the system, 
is 18. This would correspond to a planar 
system with m restricted to rrii = 0,-1,...,—17. How¬ 


ever, traditional numerical diagonalization identify the 
highly-gapped iz = 1/3 state in a slightly smaller system 
where = 16 and = 0, —1,..., —15. The numeri¬ 
cal ground state is a state with M = —45 and exhibits 
the signature of a quantum Hall state in numerical tri¬ 
als: a non-degenerate, translation and rotation invariant 
ground state with a strong energy gap. This numerical 
ground state is nearly identical to the famous Laughlin 
ansatz wave function for many different numbers of par¬ 
ticles [3, [nil and has been identified as the ground 
state of the zz = 1/3 system. 

The small correction to the filling factor calculated us¬ 
ing Eq. o is due to the finite size of the system, and 
the uncorrected filling factor approaches the ideal ratio¬ 
nal fractions of the experimental system in the thermody¬ 
namic limit. The precise locations of many quantum Hall 
states have been established in numerical trials for a wide 
variety of states. The M of the Laughlin filling functions 
(iz = 1/m for m = odd integers) are easy to establish 
based on the form of the Laughlin wave function. For a 
system with N particles, the Laughlin 1/m wave func¬ 
tion on the plane always occurs at M = —mN{N— l)/2. 
The relative azimuthal angular momentum, M, of in the 
independent particle picture is always a good quantum 
number, and is the same as the M of the hyperspherical 
picture. As a result, we use the conventional system to 
identify which M manifolds in the hyperspherical system 
contain the previously identified quantum Hall states. 

The locations of the Jain composite fermion states on 
the plane (e.g. v = 2lb,'ill,etc.. ) were established 
by using the Jain composite fermion picture ii. The 
composite fermion sequence is found for choices of v* = 
1, ±2, ±3,... and positive integer p at the filling factors 
V given by 


The strongest composite fermion states correspond to 
smaller values of |zz*| and p. We have used the com¬ 
posite fermion construction on the Haldane sphere Q to 
identify the planar M values for the Jain states. Because 
these electronic wave functions on the sphere involve only 
single-particle wave functions in the lowest Landau level, 
they can be mapped straightforwardly from the Haldane 
sphere to the infinite plane according to a stereographic 
mapping [53 |. The planar M values for the strongest 
quantum Hall states for three, four, and six particles are 
shown in Table m The filling fractions of the composite 
fermion picture {vcf) are corrected to their values in the 
thermodynamic limit. For a more general system, the 
values of M for the strongest composite fermion states 
can be calculated according to Table IHII Other hyper¬ 
spherical filling factors cannot be matched to a filling 
factor in the thermodynamic limit in the absence of ei¬ 
ther a theoretical picture (i.e. CF theory) or a series 
of numerical trials with many more particles that would 
allow the unidentified states to be extrapolated to be ex¬ 
trapolated to the many-particle case. 



















